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A self consistent effective field theory of modified gravity has recently been proposed with spon- 
taneous breaking of local Lorentz invariance. The symmetry is broken by a vector field with the 
wrong-sign mass term and it has been shown to have additional graviton modes and modified dis- 
persion relations. In this paper we study the evolution of a homogeneous and isotropic universe in 
the presence of such a vector field with a minimum lying along the time-like direction. A plethora 
of different regimes is identified, such as accelerated expansion, loitering, collapse and tracking. 



I. INTRODUCTION 

The evolution and current state of the universe has 
become amenable to a number of astronomical and ex- 
perimental probes. In the past few years, it has be- 
come possible to measure its geometry, expansion rate 
and constituents with unprecedented precision. The re- 
cent measurement of the anisotropy and polarization of 
the cosmic microwave background with the Wilkinson 
Anisotropy Microwave Probe satellite has further im- 
proved constraints on cosmological parameters We 
now believe that we have identified an accurate and con- 
sistent model of the universe based on the general theory 
of relativity. 

However our current model relies on the existence of 
extremely exotic forms of energy: dark matter which 
clumps gravitationally but does not interact with light 
and dark energy which is repulsive under gravity. Fur- 
thermore, we must posit that 95% of the energy density 
of the universe is taken up by these exotic components. 
There are proposals for their fundamental origin but as 
yet no compelling explanation. 

An alternative possibility is that gravity is not what 
it seems and that general relativity should be modified 
on certain scales. A notable example of this is the rel- 
ativistic theory of Modified Newtonian Dynamics 0, 13 • 
In this theory, dark matter is replaced by a modification 
of the gravitational interaction, through the inclusion of 
a scalar and vector field. These extra fields can com- 
pensate for the absence of cold dark matter on galactic 
and galactic-cluster scales, and at the cosmological hori- 
zon. Other proposals involving additional metrics and a 
combination of scalar, vector and tensor fields have been 
proposed. Typically these models are constructed with 
the goal of reproducing observation but with little or no 
basis on fundamental principles. Mostly they are plagued 
with inconsistencies at the quantum level (see for exam- 
pie i). 

Recenty, one of us has proposed an action for Ein- 
stein gravity coupled with a vector field that leads to 



modifications to conventional gravitational interactions 
P|. The action is constructed strictly according to the 
rules for an effective field theory with a mass scale cut- 
off, M: all Lorentz invariant terms containing the vec- 
tor field and metric up to a predefined order in E/M 
(where E is the energy scale at which the theory is be- 
ing considered) are included. The theory is manifestly 
self-consistent at energy scales below the cutoff M, and 
the usual problems of theories of modified gravity, such 
as ghosts 0, strong-coupling 0, and discontinuities @ 
are absent. The vector field has a Lorentz-invariant mass 
term with the wrong sign, which leads to a vacuum expec- 
tation value for the vector field and spontaneous breaking 
of Lorentz invariance. As a result one finds that there are 
additional graviton modes and modified dispersion rela- 
tions. Moreover, the new terms in the action lead to 
additional terms in the equations of motion which affect 
the dynamics of the universe on very large scales. It is 
this latter aspect of the theory that we wish to explore 
in this paper. 

Our work in this paper complements that presented in 
where it has been shown that Lorentz- violating vector 
fields may slow down the expansion rate of the Universe. 
In their analysis, the authors consider a fixed norm vector 
field, pinned down through a Lagrange multiplier term in 
the action. Furthermore, t hey examine a simplification 
of the theory presented in [l(| by restricting themselves 
to quadratic terms in the the vector field. In this paper 
we include higher order terms and allow the norm of the 
vector field to vary. A fixed norm case will be consid- 
ered as a limiting case where a coupling constant in the 
potential for the vector field becomes infinite. 

The layout of this paper is as follows. In Section UTI we 
present the action and briefly discuss other attempts at 
studying Lorentz violation from vector fields. In Section 
IIIII we write down the general equations of motion. We 
then specialise to the case of a homogeneous and isotropic 
spacetime. As will become apparent, the parameter space 
of this theory is immense. There are effectively 6 coupling 
constants and three energy scales that play a role. We 
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briefly discuss the range of parameters that need to be 
explored in I1VI In the following three sections we look 
at the dynamics in more detail. In Section E] we focus 
on the case in which the norm of the vector field is fixed. 
This corresponds to the case where a coupling constant 
of the potential is infinite and the theory becomes an 
extension of the theory studied in 0, ^| . In Section [22 
we set the potential energy to zero and allow the vector 
field to roam freely. This allows us to study the impact 
of the multiple derivative couplings between the vector 
field and the metric. We then look at the more general 
case in Section lYlII and conclude in Section lYlIII 



II. AN EFFECTIVE THEORY OF 
SPONTANEOUSLY-BROKEN LORENTZ 
INVARIANCE: A RECAP. 

The effective action for gravity coupled to a vector field 
can be written as 



S = S G [g,A] + S m , 



(1) 
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+c 3 S a n S b m + c 4 A a A b g mn , 
K(2) ab mn = c 5 S a n A b A m + c 6 g ab A m A n 

+c 7 S a m A b A n + c 8 A a A b A m A n 

and the dimensionful coupling constants, Cj are related 
to the dimcnsionlcss ones through 
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As compared to the theories studied in 0, 0] we are 
not restricting the vector field to have fixed norm, nor 
are we restricting the vector field to be exactly timelike. 



III. COSMOLOGICAL EVOLUTION 
EQUATIONS 

We can find the equations of motion for the effective 
theory in a general form using the formulation ofJ3}. The 
Einstein equations are 



Gab — T ab + 8nGT^l, 



(5) 



The gauge coupling constant is small (j € 1) and the 
dimensionless coupling constants, {a\,a%, /?i,/?2,/?3,/34,7) 
and n a n a are of order unity. Note that n a n a is to be fixed 
during variation. The ellipsis denotes terms suppressed 
by powers of M or Mp^; note that such a theory leads 
to corrections to General relativity of order M 2 /Mp L in 
the weak- field limit. The S m contains all other fields that 
contribute to the total action. 

We can rewrite the action Sq in a form akin to that 
presented in as 

S G [g,A] = ^ J d 4 x(-g) 1 *[R + K ab rnn V a A m V b A n 
—j8irG(A a A a — M 2 n a n a ) 2 ], (4) 
where the kinetic kernel is defined through 



K(l 



> ab 



K(ir b mn + K(2r b m 

cig ab g mn + c 2 S a m S b n 



where G ab is the Einstein tensor and T" b l is the energy- 
momentum tensor of the fields contained in S m .T ab is 
defined to be 
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in which 



jb _ V-ab tj An 
1 n ~ ly ran v a^ 1 



7 i h 



inn 



K = K ab mn V a A m V b A n , 



'V = (g ab A a A b - M"n a n a ). (7) 

The vector field equation of motion is given by 

= VJ a m + V a J a m - 2 Ci A e V e A a V m A a 

-c 5 S a n A e V a A e V m A n - c 5 S a n A b g ern V a A m V b A n 
-2c 6 g ab A n g em V a A e V b A n 
-c 7 S a e A n V a A e V m A n - c 7 g nm S a e A b V a A e V b A n 
-2c 8 A b A a A n V m A a V b A n 
-2c 8 g em A a A b A n V a A e V b A n 

+32rrG 1 VVg ma A a . (8) 

We now wish to restrict ourselves to a homogeneous 
and isotropic universe, i.e. one in which the metric is of 
the form 



g = -dt 2 + aft^ydx' <g> dx J 



(9) 



(where we have restricted ourselves to Euclidean spatial 
sections). To do so we must pick n a n a = —1, i.e. the 
minimum of the vector field must be time-like(A spacelike 
minimum would require the vector field to have spatial 
components and would violate isotropy) . The vector field 
then has the isotropic and homogeneous form 

A = (A(i),0,0,0). (10) 

The contribution from S m is taken to be: 

T m = pU<8>U + P(g + U<8>U), (11) 

where g a bU a U b = -1. 

The equations of motion now simplify dramatically. 
The non- vanishing Christoffel symbols are 

1 11 — 1 22 — 1 33 — aa i 

T^l T^l T^2 T^2 T^3 T^3 

1 01 — 1 10 — 1 02 — 1 20 — 1 30 — 1 03 — ~J 

which can be used to calculate various contributions to 
equations and |S1 for example 



V e A e = A + 3- 
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iA 2 -M 2 n°n°), 



K = {c 1 +c 3 ){A 2 +i-A 2 ) + C2 (V e A") 2 
or 

-(c 4 + c 5 + c 6 - c 8 A 2 )A 2 A 2 
-c 7 A 2 Av e A e , 

and the dot denotes differentiation with respect to t. 



The Friedmann equation is now modified. From the 
00th component of the Einstein equations we have 



= 8nGp + ^ c,A 4 
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(12) 
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From the trace of the spatial part of the Einstein equa- 
tions we find 



2 j = 87rGP + VciTi 

-4ttG7F, 



where 
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The equation of motion for the vector field becomes 

= Yl Ci ** ~ IfoG-fWA, (13) 
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where 



*i = A + 3-A~3^A, 
a ar 

* 2 = 3-A~3^A + 3-A + A, 



a ■ Or . 
A + 3-A-3—A, 
a a 1 



a 



4 = -3- A 2 A - AA 2 - A Z A, 
a 



A 2 



* 5 = —3-A A - AA - A A, 
a 

* 6 = -3-A 2 A- AA 2 - A 2 A, 
a 

^ 7 = +l-A 3 + 3A 3 ^ + 3-A 2 A-AA 2 -A 2 A, 



2 a 



s = 3-A 4 A + AA 3 A 2 
a 



A A. 



The additional terms significantly complicate the equa- 
tions of motion. In particular, we now find that the 00 
equation (a constraint equation in general relativity) con- 
tains second derivatives of the scale factor which appear 
in T ao . This situation is similar to the situation one en- 
counters in higher-derivative theories of gravity. It is not 
unexpected: in Jj|], two additional graviton modes are 
found, a clear indication that we should find more degrees 
of freedom in this theory than are present in conventional 
Einstein gravity. 



There are also three relevant dimensionful scales, 
namely the Planck scale Mp, the cutoff M, and the ob- 
servation scale Mq. From these, we define the two dimen- 
sionless ratios r = M/Mp and tq = Mo/Mp. Corrections 
to General Relativity in the weak field limit are of order 
r 2 . The observation scale can be defined in terms of the 
energy density, p by fixing the scale factor, a to be unity 
at the time of observation. We then have 

a n 

where n depends on the equation of state of the mat- 
ter component. So, in addition to the six dimensionless 
parameters we have two dimensionless energy scales. 

We clearly have a very large space of parameters to ex- 
plore. To do so, and with the aim of illustrating qualita- 
tively the possible cosmological dynamics, we will restrict 
ourselves to considering subspaces of the full parameter 
space in what follows. To this end, in the rest of this pa- 
per, we will try to identify what kind of behaviour each 
of the terms in the effective action, pinned to a given 
parameter, will generate. 



V. STRONG COUPLING LIMIT: 7 -» 00 

In the limit of 7 tending to infinity, we expect the vec- 
tor field to be fixed at the minimum of the potential. This 
can also be realized by including a fixed-norm constraint 
to the action, of the form 
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IV. THE PARAMETER SPACE 

Because of the large number of independent, Lorentz- 
invariant, two-derivative terms in the effective action, the 
parameter space for this theory is large. Even though we 
have restricted ourselves to a homogeneous and isotropic 
universe, we must still contend with eight coefficients, Cj, 
and the coupling constant, 7. Inspecting the equations 
of motion, we find that there are degeneracies. C\ and C3 
multiply equivalent terms as do C4, C5 and Cq. Putting 
them all together and reverting back to dimensionless 
parameters we find that we are left with six independent 
parameters. They can be organized into various groups 

Kinetic terms, quadratic in A: ot\ and a 2 

Kinetic terms, quartic in A: fa and fa 

Kinetic terms, sextic in A: fa 

Potential energy terms: 7 

Unless an additional symmetry principle is proposed, 
none of these terms can be discarded. In order to illus- 
trate the possible dynamics, we shall restrict ourselves to 
sub-spaces of the full parameter space in what follows; 
the full dynamics are simply too complex. But we em- 
phasize that the true spirit of effective field theories does 
not allow us to selectively discard terms. 



16ttG 



d 4 x(-g)2\(A a A a ~ M 2 n a n a ), (14) 



where A is a Lagrange multiplier, which yields the con- 
straint upon variation. Now, the time derivatives of A 
vanish. If n a n a = — 1, as if the norm of a timelike unit 
vector, then we have A=(M, 0,0,0). The equations ob- 
tained by varying the action with respect to the vector 
field and the 00th component of the metric then reduce 
to 

= 3^( Cl +c 2 + c 3 + c 7 M 2 ) (15) 
-3^(c 2 -ic 7 M 2 )-A, 
_ g 7TGp ^^M 2 (c 1 -c 2 + c 3 + 2c 7 M 2 )(16) 



+3 ^M 2 (c 2 -i C7 M 2 ) 



Eliminating the Lagrange multiplier A and converting to 
a and (3 coefficients yields the following modified Fried- 
mann equation 



+ 2 _ ai ) r 2 )3 « =87rGp . 
2 a z 



(17) 



In this limit, we then have an effective rescaling of the 
gravitational constant G in the cosmological background, 
as found in |9(. 
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VI. WEAK COUPLING LIMIT: 7^0 

Let us now look at the dynamics of the system which 
result from discarding the potential term. The vector 
field is now free to vary subject to the couplings between 
its kinetic terms and the metric. 

We first focus on the case in which a\ 7^ 0, with the 
remaining coupling constants vanishing. Varying with 
respect to the vector and wth component of the metric 
yields 

3A- = 0, 

a 

'a /d x 
2-+ - 
a \a 



8irGP + a^nGlAA + 2A 2 - 



+ (£) A 2 + A-AA + A 2 ]. 
\a J a 

The first of these equations has the solution 
a(t) = —(t + ti— — ), 

*i Si 

where at and Sj are the intial conditions at time U. Note 
that we can always rescale t — ► t+U 1 -^, such that acct, 
and we will do so from now on. 

To solve the second equation we define a new variable, 
X = \a 4 A 2 . We then have 

X = 2a 3 aA 2 +a A AA, 

X = a i AA + 2a 3 hA 2 + a 4 A 2 + 6a 2 a 2 A 2 + 8a 3 aiA 

We can now rewrite the second equation as 
2 



X-^X + 6 [ a - 



X = 



I 



{a 2 a 2 + 8irGPa 4 ), 



whose homogeneous part has the solution 

X = At 2 +Bt 3 . 

The dominant term at late times comes from the partic- 
ular integral 



X 



WnGaitf 



The particular integrals due to the pressure term of the 
matter/radiation are 

X = — ^-< 2 lni, radiation era, 
5ai 

X = — t lni, matter era. 
ai 

We can find the vector field from the definition of X. The 
dominant solution is 



As we can see, in this reduced space of parameters we 
have an attractor solution given by a oc t and A = 

y/l/2 ai Mpi. 

An altogether different type of behaviour emerges if we 
consider all constants to be zero, except for oli- The two 
equations to solve now are: 



A + 3-A-3(-)M + 3-A = 0, 
a a a 

2- + = -8nGP-8irGa 2 {AA + 3A 2 ^ + -^) 2 A 2 

a \a ) a 2 a 

+6-AA+±A 2 ]. 
a 2 

Defining Y = a 3 A,we find that the first equation reduces 
to 



Y-3-Y = 0, 

a 

with solution 

Y = Ca 3 , 
Y = C I dta 3 + D. 



We can also rewrite the second equation in terms of Y, 
which conveniently simplifies to 



2-+(-) = -8nGP - A-nGai^r 
a \a I a° 



Let us consider the case where P is negligible. Then, 
with b = In a, we find 

b + h 2 = -2nGa 2 C 2 = J, 

which can be solved as follows. Firstly, we consider sepa- 
rately the cases where J is positive or negative, denoted 
by J+ and J_ respectively. Integrating, we find 



and 



h= -J^Ltan(i/-^-(t-t )), 



respectively, where t a is a constant of integration. Inte- 
grating again yields 



b=hn(cosh(^^(t-t a )))+D 



and 



A 2 ~ —M 2 Pl . 

2ai Pl 



b=hn(cos( ] /^-(t-t a ))) + D, (18) 
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where D is an integration constant. We now have a com- 
plete solution for a = exp(6). Note that the solution is 
singular, with a — > in a finite time, for J < 0. 

This type of behaviour arises because of the novel type 
of coupling that arises when considering a vector field. 
The Lorentz structure of the vector field leads to a cou- 
pling with derivatives of the metric and hence to sec- 
ond derivatives in the would-be constraint equations and 
equations of motion. These lead to instabilities in the 
solutions of the equations of motion, a finite time singu- 
larity in this case. These instabilities presumably signal 
deeper pathologies in the effective field theory descrip- 
tion. 

If only 03 is non-zero, the 00-th and ii-th equations 
are, respectively, 



a 



8nGp 



M 2 



a 
a 



a 

8ttGP 



A 2 A 2 - A 3 A), 



The vector field equation is 



Q 



3-A 6 A + 3A*A 2 
a 



2 

8ttG/3 3 
M 2 



A 6 A. 



2 A2\ 



-A' A 



Multiplying the w-th equation by three, and adding to 
the 00-th equation we recover 

-6- = 8irG(3P + p). 
a 

where we have used the vector equation. We see that, 
in this case, the evolution is identical to that of General 
Relativity. For the case of only /Jj ^ 0, the same applies, 
with the caveat that, in the weak coupling limit, there 
is nothing to stop A(t) approaching values beyond the 
applicability of the effective action. 

Finally, we remark that the case in which /?2 ^ is 
significantly more complicated, given the appearance of 
a wide variety of terms in each equation. We may, how- 
ever, combine the vector and ii-th equations so as to 
to obtain two equations, each containing just one type of 
second-order time derivative term. In this way, we obtain 
expressions for ^ and A. It turns out that the A equa- 
tion is singular at A = and A = (8/3ft) 1/4 \/MMp. 
We will see later that in the corresponding case with 
a natural value of 7, physically viable solutions are ob- 
tained only for positive 0%. Moreover, for acceptable val- 
ues of A the second singularity will be real and lie above 
A = M. Solving the equations numerically, we find that 
for ^4(0) < the system generically reaches the singu- 
larity at A = (see Figure [TJ, whilst for A(0) > the 
system reaches the second singularity. 



VII. THE GENERAL CASE 

We now consider solutions for natural values of the 
parameters «i,/3i,7 (order unity). Furthermore, 7 is as- 
sumed to be positive. In it was argued that r could be 




FIG. 1: Weak coupling limit response when $2 ^ and 
^4(0) < (note that the units are arbitrary). 



as high as ~ 10~ 2 and as low as ~ 10 -31 and we will opt 
to look for behaviour between these bounds. Even with 
these considerations, the parameter space remains vast 
so we will proceed, as before, by considering the effect of 
isolated non-vanishing 014 and 0i parameters. Following 
this, we will consider a few non-vanishing combinations 
of oii and /3t which yield novel behaviour. 
Reasonable initial conditions are 



6(0) 

6(0) 
A(0) 
1-4(0)1 



0, 

rM 

M, 
rM 2 . 



We will consider the universe to be radiation dominated. 

Let us first consider u\ 7^ and all remaining constants 
equal to zero. Recall that this is the case where only 
the R a }jA a A h term contributes in the action. From © 
see that, peculiarly, second derivatives of the vector field 
A are absent from the vector field equation and present 
in the metric field equation whereas 6 appears in both. 
Hence the vector equation becomes an evolution equation 
for 6 

We find for all values of r that the field A undergoes 
damped oscillations, irrespective of the sign of a±. In- 
deed, for r << 1, A is dominated by 

, = AjM(A 2 -M 2 ) 

3a 2 r 2 A + '"' 

where the ellipsis denotes damping terms, which are small 
independently of the sign of a±. In the limit of A — > M, 
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we have that b — > 0. Thus we generically approach a 
loitering solution where the scale factor a grows linearly 
in time. 

If we now consider only 7^ 0, we note that sec- 
ond derivatives of b and A appear in each field equa- 
tion, an indication of how the coupling to the vector field 
modifies the kinetic terms of the metric. Consequently, 
for instance, terms in the metric equation of order 1 or 
lower in time derivatives contribute to A to a degree sup- 
pressed by r 2 . As in the weak coupling limit [Bj], the 
sign of «2 has an impact on the evolution. It is found 
that positive values of a.2 generically lead to unbounded 
growth in A, whilst negative values (of order unity) lead 
to oscillations about M, damped roughly on a timescale 
(rM)" 1 . Therefore, at times much greater than this, the 
vector field will lie fixed at A — M, and we can read 
off the resulting contribution to the metric equation as 
Tu = (—3a2r 2 /2)a~ 2 Gu. The evolution equation for the 
scale factor now becomes: 



a a 



8ttG 



1 



|a 2 r 2 



P. 



(19) 



We can interpret this as a rescaling of the gravitational 
constant G. As discussed in Q, this can only be detected 
by comparison with regimes where the vector field may 
rescale G by a differing amount. 

It is interesting to note that another consistent solution 
may emerge for high values of r- that is values of order 
unity. In this range we see that A may be attracted to a 

fixed solution at A = . / _g g Mpi = B, now close to M 



(see top panel of Figure |2J - At A — B we have that 

8 = ~4M^ (3B2 " M2)[B2 ~ M2) ' 
= C. 

To understand the nature of this solution, we consider the 
evolution of the ratio of the scale factor a(t) to the scale 
factor aGn(t) resulting from the same initial conditions 
but in the absence of the vector field, i.e. as in pure 
general relativity. This ratio is plotted in the bottom 
panel of Figure [3 The stability of this solution depends 
on the sign of C: For C > 0, b grows linearly with time 
and A is constant, leading to an asymptotic solution at 
late times of the form 



o(t) 
A(t) = B, 



a c{t-t a f 



(20) 



where t a is a constant of integration. In this limit the evo- 
lution of the vector field and metric are overwhelmingly 
dominated by terms in b. We can consider the stability 
of this solution by considering the evolution of small per- 
turbations to b and A. The first order perturbation A\ 
to A is found to obey 



A 1 +3CtA 1 +9C 2 t 2 A 1 = 0, 
The general solution to this equation is: 

„ . . , e^- Tlr ,n/3 1 iVSCt 2 ^ 
A x {t) - hx—^-Wj^-^-,-,-^—), 



(21) 



+k 2 — W w (^ r ,-,^^), (22) 
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where ki and &2 are integration constants and Wm and 
Ww are Whittaker M and W functions. We have checked 
that A± — > as t —* oo for the appropriate values of r, 
such that the solution is stable. 

In the case that C < 0, b will decrease linearly in time 
whilst the vector field is in the vicinity of B. Note that 
with the opposite sign of C, perturbations (as in J2J) are 
no longer stable. Thus we may expect that if the system 
can reach a fixed solution at B it will only do so briefly 
before moving away. This is amply illustrated in figure 
13 For r = 1 (C > 0) the system quickly settles to the 
solution A = B prompting runaway growth in cl/olcr- 
For r = 0.8 (C < 0) the system settles briefly at A = B, 
during which time a/acR is enhanced significantly before 
departing to the solution A = M. For r = 0.6 (and all 
lower values of r) the vector field never reaches B and 
settles to the A = M tracking solution described by (|19|l . 

We now restrict ourselves to 02 7^ and discard all 
other terms. We note that for the ranges of r considered, 
the values at which A is singular (see I VI (I arguably lie 
beyond the regime of validity of the effective action. For 
negative /?2, A experiences unbounded growth in the di- 
rection of the initial perturbation, as in Section IVII and 
so will generically reach singularities. For positive 02, 
the vector field undergoes oscillations about M. Numer- 
ical exploration suggests that the system will evolve so 
that the oscillations will be of increasing frequency and 
amplitude up until t « (Mr) . Terms in T then act 
to halt the expansion of the universe leading to eventual 
collapse (see Figure EJ. 

If we consider /3 3 ^ or f3± ^ 0, we recover familiar 
behaviour. The vector field's stress energy tensor con- 
tributes only first derivatives in time and terms such as 

■ 2 

^ and ^ do not appear in the vector field equation; thus, 
A contributes and evolves much like a scalar field. It is 
found that for positive values of 03 and for negative val- 
ues of /?4, the vector field oscillates about A. It undergoes 
Hubble friction, the oscillations decaying on a timescale 
~ M _1 . Settling towards A = M, the potential terms 
vanish as before and so the effect on the expansion of the 
universe is as in the weak-coupling limit i.e. identical to 
that of general relativity, 
w 

An interesting simplification is obtained for particular 
the combination of coefficients, given by ot\ — &2 < 0. 
With this combination, b drops out from the vector equa- 
tion. For acceptable values of r, the vector field oscillates 
about M, the background expansion gradually settling 
the field to this value. This happens on a timescale of 
roughly M . Thus at times larger than this the metric 
equation reduces to 



a a 1 — 2r 2 a\ 

Again this behaviour can be interpreted as a rescaling of 
G or, in this the tensor T tracking any form of 

matter field in the universe. 




0.8 - 

a o.6 \ /Vv»» 

- 

0.4 - 
0.2 - 
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t 

FIG. 4: Comparative evolution of A for r values in ration 
10:6:1 (solid:dashed:dotted) with r = 1 being the highest 
scale. Such a high value is chosen only for ease of illustra- 
tion. ^ = —1,7 = 1, and M p i=l. 

Finally, we would like to make a general point about 
the appearance of singularities in this system. The field 
equations can be written schematically as: 

i = i(...), 

where / is of the form 

/ = (-2 + r 2 (2 ai -3a 2 )A 2 ) (23) 
x(~a 2 + (P 3 + P2)A 2 -(3 i A 4 ) 

-{ai - a 2 H —J 

x(3a 1 -3a 2 -^ 2 -A 2 )r 2 ^ 2 , 

and where A = j^. When the terms in parentheses are 
nonvanishing, there will may occur singularities in the 
evolution for particular values of A . An example of this 
was encountered in IVII for only @2 ^ 0. 

We expect r to be of order 10~ 2 or smaller and A 
to be of order unity. Hence, generally, (—2 + r 2 (2ai — 
3a2)A 2 ) ~ —2. Therefore, the second terms in paren- 
thesis may be expected to be suppressed relative to the 
first by r 2 . It is unlikely then that / may vanish due 
to equality of the first and second terms in parenthe- 
sis. An alternative is for both groups of terms to vanish 
identically. This would seem to require of the first that 
8 = — «2 + (/?3 + /?2 ) ^4 2 — 04 A 4 vanishes along with at least 
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one term from the second group. When considering iso- 
lated non-vanishing coefficients, restrictions on their sign 
were found to be a2, ft < and /?2,/?3 > respectively. 
This would seem to imply that S is inherently positive. It 
may only vanish when each of the coefficients are zero, in 
which case / may only vanish if a\ is alzo zero. However, 
we emphasize that the restrictions on the sign of individ- 
ual coefficients need not hold when general combinations 
of «j, Pi are considered. 

VIII. CONCLUSIONS 

In this paper, we have studied the cosmology of the 
model proposed in || . As expected, there is a wide range 
of possible behaviours. We achieve accelerated expansion 
if we have a term with a 2 ^ and a high mass scale M. 
The system will scale when a± = 0% or when a 2 ^ with 
a low mass scale. If /?2 ^ 0, the universe will recollapse, 
while it will loiter if we only have ct\ ^ 0. Clearly we 
have only looked at a small subset of the parameter space 
but our analysis has allowed us to probe these different 
regimes. It also allows us to make some comments about 
the structure of the theory. 

Firstly, we stress again that this is an effective field the- 
ory, valid at energies below the cutoff, M. In some cases, 
we have found that this theory is unstable, in the sense 
that runaway solutions push the theory beyond its regime 
of validity. To properly understand its behaviour in these 
regimes, we would need to at least consider higher order 
corrections, and ultimately some ultra-violet completion 



of the theory would be necessary. Such a completion 
might conceivably simply involve extra fields or, 
more likely, more complicated dynamics. 

Secondly, the Lorentz structure of the vector field in 
this framework introduces a novel phenomenon. For the- 
ories with C2 7^ and cj ^ we find second derivatives of 
the scale factor appear in what would have been (in stan- 
dard General Relativity) constraint equations, and in the 
evolution equations for the vector field. We find that this 
leads to possible instabilities in the cosmology. Clearly, 
there is a need for a complete perturbative analysis of 
this theory in a cosmological setting, akin the study of 
Gauss-Bonnet or higher-derivative modifications to grav- 
ity [12. 

Thirdly, the results found here can be added to known 
cosmological consequences of Lorentz violating fields: de 
Sitter expansion and dust-like stress energy tensor in 
background 01 , pervasive tracking solutions in back- 
ground , and generating the instability permitting the 
growth of large scale structure 0] . 

Finally, the theory discussed here, though complex in 
its structure, should be amenable to a detailed compar- 
ison with current cosmological observations. We have 
laid down the framework for looking at constraints on 
the background evolution. The next step is to construct 
the evolution equations for linear perturbations. With 
these in hand, it should be possible to harness the wealth 
of new, high precision cosmological observations and use 
them to provide constraints on this theory. 
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